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The Longuet-Higgins-Berry's phase has remarkable consequences for charge transport in molecular 
rings. For generic (conical) crossing, where the phase is n, a vanishing cause can lead to a diverging 
response in the amount of charge transport. Away from level crossings, when the phase is 0, a 
vanishing cause leads to a vanishing response. The divergence of the response near crossing is related 
to, but distinct from, the divergence that occurs in the generalized susceptibility. We illustrate this 
behavior for quantum models of molecular rings driven by a running wave of small amplitude at 
zero and finite temperatures. 
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Consider a molecular ring, such as a benzene or a 
triangular molecule X3, such as Na3. In the Born- 
Oppenheimer limit of heavy nuclei one can consider a cy- 
cle of deformations where each nucleus is displaced only 
slightly from its initial position and eventually returns to 
it. The question that we want to focus on is: What is 
the electronic charge transported around the molecular 
ring in one such cycle? As we shall explain below, there 
are two cases: If the cycle of atomic deformations can be 
shrunk to zero without trapping a point of level cross- 
ing (of the electronic energy levels) then one gets normal 
behavior in the sense that the weaker the deformation, 
the less the charge transported in one cycle. If, however, 
the cycle of atomic deformations pinches a point of level 
crossing, then the smaller the cycle, the larger the trans- 
ported charge. We dub such anomalous behavior, where 
the weaker the cause the larger the effect, homeopathic. 
We shall illustrate this for Hiickel (tight binding) models 
of molecular rings at zero temperature. 

Homeopathic charge transport is intimately related to 
the Longuet-Higgins Q and Berry's phase M: For time 
reversal invariant Hamiltonians level crossing is reflected 
in the sign of the electronic wave function undergoing 
a cycle of deformation. This sign is —1 if the cycle 
pinches a (generic) crossing and 1 if it does not. The 
Longuet-Higgins phase has important consequences for 
(molecular) rovibronic spectra in Born-Oppenheimer and 
Jahn- Teller theory and plays a role in molecular dy- 
namics ||. The theory of adiabatic transport adds the 
observation that Longuet-Higgins phase has direct con- 
sequences also for electronic properties, and not only to 
molecular properties. In particular, a phase —1 implies 
homeopathic charge transport. There is both theoretical 
and experimental evidence that molecular trimers such 
as Na3 Q and sym-triazine as well as other trimers 
and systems have —1 Longuet-Higgins phase. 

Adiabatic charge transport near gap closures for infi- 
nite chains has been studied in For these there is 
no homeopathic divergence. Instead, a vanishing cause 
can lead to a finite and quantized charge transport when 
the Longuet-Higgins phase is —1. The charge transport 
in finite molecular rings is, therefore, more singular than 



that in infinite chains. 

The homeopathic behavior that occurs for (out of equi- 
librium) charge transport is related to, but distinct from, 
the divergences that can occur in thermodynamic equi- 
librium of generalized susceptibilities at T — 0. A nec- 
essary condition for either, at least for the simple model 
systems we consider below, is that quantum energy levels 
cross. But, while thermodynamic susceptibilities probe 
the singularity of the energy surface near crossing, non- 
equilibrium adiabatic transport probes the singularity of 
the surface of eigenstates near crossings. It is possible 
for one of these surfaces to be singular without the other 
being singular. We shall return to this issue below. 

As we shall see, finite temperature introduces a cutoff 
of the homeopathic divergence. In some cases, a Jahn- 
Teller instability can censor the homeopathic divergence 
even at T = 0. 

Consider, for simplicity, the Hiickel (tight-binding) 
model Hamiltonian for non interacting electrons in a gen- 
eral triangular molecule of three identical atoms. Al- 
though the example of a molecular trimer is special, it 
turns out that it describes the generic situation near level 
crossing. The Hamiltonian is the 3x3 Hermitian matrix 



H(a, b, c, 4>) = E Q 




(1) 



where Eq fixes the energy scale and a, b, c are (dimcn- 
sionless, real) hopping amplitudes associated to the three 
bonds of the triangle. (The triangle is not necessarily 
equilateral.) We assume that a,b,c are all positive. We 
can, and shall, use a, 6, and c as local coordinates in the 
space of (internal) configurations of a trimer. a, 6, and 
c are actually not good coordinates globally. For good 
global coordinates see ||. However, global subtleties 
need not concern us here since we consider only small 
deformations. 

£ = exp i(f>, with <f> an auxiliary phase variable associ- 
ated with a fictitious Aharonov-Bohm flux tube which 
carries flux <f> and threads the molecule. The explicit 
form of Eq. (|l|) involves a choice of gauge for the flux 
tube. (We shall consider observables that are indepen- 
dent of this choice.) The role of <f> will become clear 



1 



below. Deformations of the molecule change the hopping 
amplitudes, and a closed cycle of deformations is a closed 
path in the three dimensional space whose points are the 
hopping amplitudes (a, b, c) G R + . Such a closed path is 
shown in Fig. 1. For notational convenience we denote 
by X the triplet (a, b, c). 

The reason for introducing the fictitious flux tube 
is to define the current operator which circulates in such 
a ring. For the choice of gauge we have made for the 
Aharonov-Bohm flux, the current is associated with a 
single bond, the c bond, and is 

/o o -*A 

(d tf)(A,0) = c£ 0. (2) 
Wo / 

This is the sole role of and in all our calculations we 
shall eventually set = 0, which is the case with no flux 
at all. In this case the Hamiltonian H(X, 0) is real and 
therefore time reversal invariant. The observable associ- 
ated to the circulating current, Eq. (^|), is pure imaginary 
when = and so odd under time reversal. Because 
the model for = is time reversal invariant, there are 
no diamagnctic (persistent) currents in any cigcnstatc. 
Let P(X,<fi) be a spectral projection for H(X,(f>), i.e. 
H(X,cj>)P(X,<t>) = E{X,cj))P{X,<t>) with E(X,(f>) S H 
an eigenvalue. The vanishing of the persistent currents 
is the statement Tr (Pd^H)\,p = o = for all X. 

We are interested in the current that flows around the 
molecule when it is slowly driven so that X traces a closed 
path in parameter space, as e.g. in Fig. 1. In the theory of 
adiabatic transport this current is related to the adiabatic 
curvature. The basic equation for the expectation value 
of the current at time t and zero flux, reads [p|-^8| : 

Tr{P t d 4> H)=Tr{Q 4>x (P)^ -X + 0{1/t 2 ). (3) 

X stands for the triplet a, b and c and the dot de- 
notes time derivative. Pt is a solution of the quan- 
tum evolution equation, with the (adiabatic) time de- 
pendent Hamiltonian H(X, 0) and with initial condi- 
tion that Pt=o is an eigenstate. P(X, 0) is an instanta- 
neous spectral projection for the instantaneous H(X, 0). 
Q<px(P) = —iP[dcf,P,dxP]P is the <j)X component of 
the adiabatic curvature r is the time scale so that 
the adiabatic limit is r — > oo. The charge transported 
around the ring in one cycle, in the adiabatic limit, is 

Q = J o dtTr(P t d^H)= j>Tr{n^x{P))dX. (4) 

Analysis of the characteristic equation of Eq. (jl|) shows 
that level crossing can occur only if a = b = c and £ = 
±1. For £ = 1, the case we consider here, the simple 
eigenvalue is 2aEo (the top state if Eq is positive) and the 
corresponding eigenvector is |0) = (1, 1, 1). The two 
fold degenerate eigenvalue is —qEq and the projection to 



its subspace is 1 — |0)(0|. For = 0, level crossings occur 
on a ray in X space. 

The essence of homeopathic behavior is the following. 
Suppose first that the cycle of deformation does not pinch 
level crossing. (For the top state this holds for any closed 
cycle in the positive quadrant.) The adiabatic curvature 
is smooth and bounded along the path and the right hand 
side of Eq. (Q) is of the order of the area of the cycle of 
deformations (by Stokes formula). Q goes to zero when 
the cycle shrinks to zero. Now suppose that the cycle of 
deformation pinches the ray of level crossing, e.g. 

a(x) = 1 + lux + uix; b(x) = 1 + x + x; 

c(x) = 1 + LUX + LUX, (5) 

where to is here a (complex) cube root of unity and x S C 
runs on a small circle in the complex plane surrounding 
the origin, see Fig. 1. As the cycle shrinks, the adia- 
batic curvature in Eq. (Q) diverges quadratically, like a 
monopole p|Jl9|], while the length of the cycle shrinks 
only linearly. It follows that now Q diverges as the cycle 
shrinks to zero, i.e. Q — O f tjt) • 

In the simple case at hand the leading divergence char- 
acterizing homeopathic behavior can be calculated ex- 
plicitly. The degeneracy splits in first order of perturba- 
tion theory, both in and in x, and the local behavior 
near crossing of the two bottom states is given by the 
2x2 matrix 

Eg f -\/30 6uix \ , , 

3 ^ Qujx V30 J ' [ ' 

This matrix has the form of a Berry spin 1/2 model and 
the adiabatic curvature of its two states is explicitly com- 
putable. So, to leading order, the equation for the charge 
Q, Eq. (||), reduces to 

Q = ±^f j A d Angle, (7) 

where Angle is the angle swept by x as it moves around 
the origin in the complex plane. The ± signs refer to 
the ground and first state respectively. A simple formula 
is obtained for a circular orbit, \x\ = const, where Q = 
■q^j- Evidently, the smaller the cycle that pinches the 
degeneracy, the more charge it transports and Q — > ±oo 
as — > 0. 

It may be worthwile to explain which aspects of Eq. (|t]) 
are general and hold for any generic two level crossing and 
what is special for the explicit model we consider. The 
overall constant ^/3/12 is special for the model. What is 
general is that the divergence scales like an inverse power 
of the distance from crossing. This can be seen be noting 
that for any (generic) two level crossing the adiabatic 
curvature diverges like the field of a monopole, i.e. like 
|ir|~ 2 , and therefore a line integral on a loop of length 
0(|x|) will scale like the potential of a monopole. 
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It should be stressed that this result does not imply 
that the current is large. In Eq. (||) only the ratio of the 
current to the rate of driving is large. The current is not 
large because, as the circle is shrunk, the rate of driving 
must also decrease in order for the adiabatic theory to 
apply. 

It may be worthwhile to point out that the charge Q in 
Eq. (Q) even though a geometric, is not a Berry's phase 
(being a line integral, rather than a surface integral, of 
the curvature). 

We now return to the question of how is homeopathic 
behavior of transport related to the divergence of the 
susceptibility near level crossing. At T = the general- 
ized susceptibility matrix x of a family of Hamiltonians, 
H(X,cj)), that depend parametrically on X and cj> is the 
symmetric matrix of second derivatives: 



Xxx(X,, 



dE 
dxidx 



■(X,<f>), xxcb(X,, 



dE 



dxi 



where E(X,(f>) is the ground state of H(X,(f>). x 1S a 
thermodynamic (equilibrium) property and as such it 
depends only on the the energy surface E(X,(f>). The 
adiabatic curvature is associated with an anti- symmetric 
matrix whose components are 



Trn xx (X, <j>) 
Trn X4> (X, <j>) 



-iTr(P[d Xi P,d Xj P]), 
-iTr(P[d Xt P,d^P}). 



(8) 



It describes transport coefficients which can not be de- 
termined from the the ground state energy (at T = 0) 
or a thermodynamic potential (for T > 0). This is evi- 
dent from the formula for the adiabatic curvature, which 
is determined by the projection P(X, </>) on the ground 
state and not on its energy E(X, (p). 

Eigenenergies and eigenstates are, of course, related. 
Away from level crossings, both E and P are smooth 
functions of the parameters. However, besides that, 
the susceptibility and the curvature are essentially in- 
dependent quantities. For example, from Eq. (Q) the en- 
ergy surface of the ground state is E((j),x) — —Eq (1 + 
v / ( /) 2 /3 + 4|a;| 2 ). So, while Q X A X ^ = °) diverges at 
the crossing, the corresponding component of the suscep- 
tibility, xx,4>(X,4> — 0) vanishes identically and, while 
Xtj>it> diverges at crossing, the corresponding curvature 
Q,^ vanishes identically. 

Finite temperature provides a cutoff to the homeo- 
pathic divergence. Two nearly crossing states, which are 
at distance e apart in energy, transport opposite charges 
Q = ±0(l/e). At finite temperature T the two nearly 
crossing states will be nearly equally populated with a 
bias of 0{e/T) towards the lower state. For the trian- 
gular molecule undergoing a cycle of deformations, the 
leading behavior of the charge transport at low temper- 
ature is: 



Q(T) = 



kEq 



where ks is Boltzmann's constant. The constants in this 
formula are special for the model. In general, for generic 
two level crossings, one can conclude that the total charge 
transport will approach a finite limit as one approaches 
the crossing so that: 



Q(T) = 0(1/T). 



(9) 



The 1/T law is reminiscent of Curie's law. 

Similar analysis can be made for a necklace of p equiv- 
alent atoms with p an arbitrary integer larger than 2. 
(For benzene p — 6.) The tight binding Hamiltonian is 
a p x p Hermitian matrix with nearest neighbors hop- 
ping only. Now suppose that such a ring is deformed 
by means of a running sinusoidal wave of commensurate 
wavelength, i.e. the hopping amplitude between the fc-th 
and (fc + l)-th atom is time dependent and is given by 
1 + 2k cos[2iv(jk/p — t/r)], where k (which is held fixed) 
stands for the amplitude of the distortion and j is a natu- 
ral number smaller than p. (This is an analog of Eq. (||)). 
When <fi = k = 0, there are twofold degeneracies at en- 
ergies 2 cos(2mn/p) for to — ±1, . . . , ±[p/2]. The new 
feature of this model is that the order of perturbation 
theory that splits the degeneracies depends on \m\, p and 
j, and can be large if p is large. As a consequence, the 
singularity of the adiabatic curvature near crossing can 
be quite strong (without the susceptibilities being singu- 
lar). The charge transported can be calculated here as 
well, but the details will be described elsewhere. The re- 
sult is that the charge transported in the pair of nearly 
crossing states near energy 2cos(2m7r/p) is: 



Q = ± 



7rq sin tylO 



pui cos(to — j /2)9 



fq-l 

n 

\k=l 



cos(m — jk)9 — cos to6 



cos(to + j/2 — kj)6 



(10) 



q > 1 is the order of perturbation theory that splits the 
degeneracy, q is the smallest natural number such that 
qj mod p equals cither 2\m\ or p— 2\m\. In the first case 
m = | to | in Eq. ( [To| ) and, in the second case, to = — \m\. 
The numerator and the denominator never vanish under 
the conditions that lead to this equation. 9 is shorthand 
for 2n/p. We see that the amplitude of the perturbation 
k enters this expression with a negative exponent: this 
is the homeopathic effect. The sign of Q reflects the fact 
that each member of the pair of nearly crossing states 
transports charge in opposite sense. Under a complete 
cycle t — > t + t the electronic eigenstate acquires a phase 
factor (-1) 9 . Eq. (|) is still obeyed. 

The molecular models discussed so far are prototype 
models of finite quantum systems. It is natural to in- 
quire what, if any, of the homeopathic behavior survives 
for infinite, macroscopic systems. A prototype of such a 
system is the infinite one dimensional chain with finite 
electron density. When the Fermi energy lies in a gap 
such a chain is nominally an insulator and the theory 
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of adiabatic transport applies. For non interacting elec- 
trons and a periodic chain, the transport behavior can 
be analyzed using standard, single electron, techniques. 
It turns out that there is no homeopathic divergence of 
transport. What survives of the homepathic behavior is 
that an arbitrarily small cycle of deformation can lead to 
a finite (nonzero) and quantized charge transport. This 
requires that the deformation pinches a gap closure. The 
quantization of transport (at T = 0) comes as it does in 
charge pumps @@ and the Hall effect §]. The theoret- 
ical framework presented above sheds light on the numer- 
ical findings of || who found divergence of the transport 
coefficients in Hubbard models of Perovskite chains. 

Until now we have only considered the electronic part 
of the Born-Oppenheimer theory. The inclusion of the 
rovibronic part is, in general, a formidable problem even 
for a molecular trimer where the intricacies of the three 
body problem come into play. We refer to || for what 
is known in general and to |^.^0| p2[ , for models. We 
shall content ourselves here with the classical limit for 
the nuclear dynamics. 

If the ground state of the undeformed molecule is de- 
generate, the total energy may be decreased by deform- 
ing the molecule. This is the classical Jahn Teller 
instability. This is the case if the elastic energy is pro- 
portional to the square of the amplitude of the defor- 
mation (i.e. is harmonic) and if the degeneracy lifts to 
first order in the deformation. If the Jahn Teller energy 
functional has a unique minimizer which breaks the de- 
generacy, then a sufficiently small cycle of deformations 
around the ground state will not encircle also the point 
of level crossings. In this case the homeopathic behav- 
ior (at T = 0) is censored at the ground state by the 
Jahn Teller instability and there is no divergence of the 
charge transport in the limit of an infinitesimal cycle of 
deformation. 
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FIG. 1. A cycle of deformation for a molecule with three 
atoms. 
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